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The paper is on the line of the path integral technique developped in preceding 
articles. Taking advantage of a semirelativistic and a full relativistic representation 
of the quark propagator in an external field we present an unified derivation of the 
semirelativistic potential and of a Bethe— Salpeter like equation for the quark- 
antiquark system. We consider three different models for the evaluation of the 
Wilson loop: the Modified Area Law model (MAL), the Stochastic Vacuum Model 
(SVM) and the Dual QCD (DQCD). We compare the corresponding potentials 
and show that they all agree at the short and the long distances. In the case of 
the Bethe-Salpeter equation we treat explicitly only the MAL model and give an 
explicit expression for the kernel. Then we show that an effective mass operator can 
be obtained which agrees with the MAL potential in the semirelativistic limit. In 
the light quark mass limit this mass operator produces straight Rcgge trajectories 
with Nambu-Goto slope in agreement with the data. Finally we briefly discuss the 
mass independence of the hyperfine splitting in the heavy-light case. 



1 Introduction 

In preceding papers0,S we have shown how, using certain path integral repre- 
sentations of the quark-antiquark gauge invariant Green function G gl (a;i, x-i', y±, 1)2) 
is possible to derive a semirelativistic potential for heavy quarks or a Bethe- 
Salpeter like (BS) equation for light and heavy quarks in terms of the Wilson 
loop integral W and its functional derivatives. Such representations are the 
consequence of corresponding representations for the quark propagator in an 
external gauge field. 

We have used two different path integral representations for the propaga- 
tor. The first one is expressed in terms of a ^ expansion and quark paths 
in ordinary phase spacm. This is appropriate for the derivation of the heavy 
quark potential. nine second one is expressed in terms of covariant space-time 
quark paths □, □, E3 and it is useful for full relativistic developments. 

Since up to now it is not possible to evaluate the Wilson loop analytically 
from first principles, to obtain explicit expressions one has to rely on mod- 
els and/or lattice simulations. In this paper we shall discuss three different 
models: the Modified Area Law model (MAL), in which ilnW^ is written as 
the sum of a perturbative term, an area term and a perimeter term EJ, 13; the 
Stochastic Vacuum Model (SVM) which is based on a cumulant expansion; 
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the dual QCD (DQCD)B, which introduces effective fields and an effective la- 
grangian which should simulate what is believed to be the real mechanism for 
the formation of the confining flux tube. For each of these models a different 
potential can be obtained. We shall compare such potentials and show that 
they agree for large and small distances, but differ in the intermediate range. 
Up today an explicit BS equation can be obtained only in the MAL case. Such 
equation is written in terms of a second order Green function H(xi, x 2 ; 2/1,2/2) 
to which G sl (x\,x 2 ;yi,y 2 ) is related for large x® — y®(j — 1,2) simply by the 
application of Dirac type differential operator. The BS kernel I(xi, j/i; X2, 2/2) 
is expressed in principle as an expansion in the strong coupling constant a s 
and the string tension a (more precisely the quantity aa 2 , a being the typical 
radius of a bound state) . In principle the method could be applied also to SVM 
and DQCD models but the actual calculations would be very complicated and 
have not been performed. In practice even the MAL kernel has been worked 
out only at the lowest order. We show that starting from the MAL BS equation 
a squared (M 2 ) or a linear (M) relativistic effective mass operator involving 
certain relativistic potential U or V can be obtained in the instantaneous limit. 
Spectrum calculation using such expressions are in progress. Obviously, in the 
scmirclativistic limit V reproduces the MAL potential. On the contrary, if one 
retains full relativistic kinematic but neglects the spin dependent terms, M 
coincides_with the hamiltonian of the Relativistic Flux Fube Model, discussed 
in Ref. £3, up to ordering prescription and in the limit of vanishing quark 
masses gives Regge trajectories with slope a' — ^7 as in Nambu-Goto string 
which are in agreement with phenomenology (a ~ 0.18GeV 2 , a' ~ 0.88). Fur- 
thermore, if the spin dependent terms are kept but the limit of heavy-light 
miarks is considered, spin symmetries similar to those discussed by Kaidalov 
Ii3 are obtained. In Sec. 2 we give the basic equation for the definition of the 
potential and BS equation; in Sees. 3 and 4 we review the three mentioned 
models for the evaluation of the Wilson loop and compare the form of the 
corresponding potentials; in Sect. 5 we report the BS kernel in the Ml 
and discuss the mass operator. The present paper is mainly based ont 
but it is in part original in the presentation and contains some new results. 

2 Definition of the semirelativistic potential and BS equation 

In the usual functional integral formulation, after integration on the fermionic 
variables, the gauge invariant quark-antiquark Green function can be written 

G^(x 1 ,x 2 ;y 1 ,y 2 ) = ^Ti(U(x 2 ,x 1 )S^(x 1 ,y 1 \A)U(yi,y2)C- 1 S^(y 2 ,x 2 \A)C). 

(1) 
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Here C is the charge coniugation matrix, U (a, b) denotes the Schwinger string 
U(b,a) = Pexpjig / dx^A^x)} (2) 

J a 

(with A ll (x) = ^\ a A^(x) and P the path ordering operator over the color 
matrices), S(x, y\ A) is the quark propagator in an external field 

(i-y»D„-m)S(x,y;A) =5 4 (x-y) (3) 

( Dfj, = — igA^(x)) and the notation 

VW) fVAM f (A)e iS ™M ' [V 

is used, Mf(A) being the fermionic determinant and SVm[/4] the pure Yang- 
Mills action. 

For a closed curve T we set across the paper (Wilson loop correlator) 



W(T) = -(TrPexp{^ dx^A^x)}) 



We are interested in loops formed by a line Ti joining y\ to x\ (the quark 
path), another line T2 joining x-i to yi (the reverse antiquark path) and two 
straight lines connecting X\ with xi and yi with yi. n ]jji term of these the 
semirelativistic quark-antiquark potential is defined b 
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at 

pdtVtx = i log W(r) - £^ jf d^^fli^ 



((D-F^x)))) - \ ± ^—T S I dx» / dx>° a] 4 



8 m j 



( ((F hl (x)F ka (x'))) - (F lfi (x)}) (F ka (x')))) , (5) 



with yi = y\ = y°,x1 — x\ = x°, T s time ordering prescription on the spin 
matrices, F^ v = ^e^paF 1 " 7 and 

U[ ))} ~ JVAe^WTiPexplig^dx^A^x)] ' U 

Eq. || is obtained solving Eq.|| for > Vo by a Foldy-Wouthysen transforma- 
tion and path integral technique and comparing the resulting expression for 



G gl with the two particle Schrodinger equation path representation. Notice 
that if SS >lv (z) = ^(dz^Sz" — dz^Sz 11 ) denotes the surface swept by the quark 
path as a consequence of an infinitesimal variation z(r) — » z(t) + 5z(t) one 
has 

^{Pexp[ lff ^ dzM M (z)]} =igP J 5S^(z)F^(z) exp{ig J dz^'A^z')}. 
From this it follows immediately 

^fe)))^-!)^^^, (i = M) (7) 
g 2 (((F^( Zl )F Xp (z 2 ))) ((F^UF^)))) = -ig ^— ^F^))). 

(8) 

etc. and so everything in the right hand side of Eq. || can be expressed in terms 
of W. Notice, however, that, both through the expression of i In W (e.g. cf. [l9]) 
and the spin dependent terms, the right hand side of ^ contains integrals on two 
separate times t\ and t 2 , while, by assumption, the left hand side contains one 
single time t. One can dispose of the above situation by integrating explicitely 
the relative time r = t\ — t% after setting 



1, ..,^1 . ... 1, 
8' 



= z 3 (t) + -(-ly+Wz^t) + -r%(t) + ... (9) 



keeping only the appropriate order term and eliminating the second time 
derivatives by partial integration (notice that z° = = 0(z 2 )). 
On the contrary, if one sets 

S(x,y;A) = (i^D^ + m)A(x,y;A), (10) 

Eq. U gives 

(D,D» + m 2 - \o» v F^) A(x, y; A) = -S 4 (x - y) . (11) 
Then, for large time intervals one can write 

G gi (xi,x 2 ;yi,y 2 ) = («7i^im + mi)(ij 2 d 2u + m 2 )H(xix 2 ; yij/a), (12) 

with 

H{x 1 x 2 ;y 1 y 2 ) = ^Ti(U{x 2 ,x 1 )A^(x 1 ,y 1 ; A)U( yi ,y 2 )A (2 \y 2 ,x 2 ; A)) (13) 



where the tilde denotes transposition of the color matrices. Furthermore, let 
us assume that ihxW(T) can be set in the general form 

i]nW=^y] f ' * d n [ ' d^Eizuz'^Zi,^) + ... , (14) 
1 i,j Jo Jo 

where the dots stand for A 2 terms involving similar fourth order integrals, A 3 
terms involving six order integrals and so on (A being a parameter that has been 
introduced for convenience) and E is an homogenous function of degree in z\ , 



Z2- Then, using 
representationQ. 



(x, y\ A) the covariant path integral Feynman-Sphwinger 
13, □, one can obtain the inhomogeneous BS equation] 



H(xi, x 2 - y lt y 2 )^H^ 1 \x 1 ~y 1 )H^(x 2 ~y 2 ) - i J d 4 ^d^ 2 d 4 m d 4 m 

H^( Xl - a) h^(x 2 - &) /(&, 6; m, m) H(r,u m; yi, 2/2) (15) 

with the kernel / expressed as expansion in A and H(x — y) = (A(x, y\ A)). At 
the first order in A one has 

J(£i,£2)»7i)»72) = -4 J d ^^ R ( — ^ ^ ??2 ; fc i ; fc 2)exp{-^[(^i-??i)fci + fe-7? 2 )fc 2 ]}, 

(16) 

with R defined by 

dn dT 2 E(z 1 ,z 2 ,p 1 ,p 2 )(So 1 Sn 2 )~ 1 = / dr\ \ dr 2 R(z 1 , z 2 ,p 1 ,p 2 ) 

(17) 

having set = T s exp [ - \ f° dra^ SSl ?„ {z) ] ■ 

3 Models for Wilson Loop Evaluation 

3. 1 Modified Area Law Model 

As mentioned, MAL model consists in assuming i In W to be the sum of its 
pcrturbative value, area and perimeter term 

i log<W(r)) = i log(W(r)) pert + aS min + x -cp , (18) 

where obviously >!? m i n denotes the minimal area enclosed by IhaJoop T and P 
its length as suggested by the Wilson strong coupling limit Ei, 0, while a and 



C are in practice treated as independent adjustable parameters (which must 
however agree with the lattice simulations). At lowest order one has 

i{\nW) pcrt = -9 2 / dh / dt 2 zZz» D^( Zl - z 2 ) . (19) 

Furthermore 

i 

' du^ du^ \ / 9u M du^ \ ( du^ du^ x 



r dt f . 

Jti Jo 



Smin — min / dt da 



dt dt J \ da da J \ dt da 

^ (20) 

P = J2 r° dt^z^ , (21) 
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where = u^(t,a) with y° < t < x°, < a < 1 is the equation of an 
arbitrary surface enclosed by T and satisfying therefore boundary conditions 

u"(t,0)=^(t), u"(t,l) = *?(*)■ (22) 

In practice an equal time straight line approximation is commonly adopted for 
-Smiii- This amounts to setting 

u°(t,a)=t, u(t,a) = a zi(t) + (1 - a) z 2 (i) , (23) 

Smin= [ dtr [ da[l-(ai 1T + (l-a) i 2T ) 2 ]i = (24) 
Jti Jo 

= r dn r 2 rfr 2 «5(z0-z 2 °)|z 1 -z 2 | f ^{z0i 2 °-[az 1T i 2 + (l- a )z 2T i0] 2 }^ 
Jo Jo Jo 

with z^ T = (S hk — f h r k )zj (obviously in the first step zj = in the second 

one Zj — -jf-)- This equation it is exact to order 1/m 2 and for particular 
geometries. The second step is given in terms of a covariant parametrization 
of the quark and the antiquark paths and is useful in the derivation of the BS 
equation; however it is not Lorentz invariant and it is assumed to be true in 
the center of mass frame. 

3.2 Stochastic Vacuum Model 

Using the non abelian Stokes theorem and cumulant expansions, one can write 
(W(T)) = (p exp(ig J s dS^(u)F^(u,x )^ (25) 



dS^ Ul (mi).. 



3=1 

dS 1 **"* (ujXF^^uifXo) ■ ■ ■ F H ^(uj,x a )) cum .(26) 

where F^(u,xo) = U(xo,u)F pl ,U(u,xo) and the cumulants (. . .) CU m are de- 
fined by 

(F(l)) cum = , (F(l)F(2)) cura = (F(1)F(2)) - (F(2)) , ... 

(27) 

5 is an arbitrary surface enclosed by T, xq a reference point on S^-. Ps an 
ordering prescription on S and the U(u,xo) are the Schwinger strings@. 

The basic approximation consists in assuming that the second cumulant 
is dominant and actually independent of xq. Then, since the first cumulant 
vanishes trivially, one can write 

\og(W(T)) = - 9 - J dS^{u) J s dS x P(v)(F^(u,x Q )F Xp (v,x Q )) cum . (28) 

and, taking into account Lorentz invariance, 

g 2 {F pu (u,x )F\ p (v 1 x )) cum = g 2 (F p „(u, x )F\ p {v, x )) 

{1 d 
{^\5 up - 5 w 5 vX )D{{u - v) 2 ) + 2^g^~^ u ~ v )^»p - ( u ~ v) p 5 v \) + 

+ ^-((« - u )p<W -( u ~ v)xS PP )]D l {{u - «) 2 )| (29) 

where (3 = fg ^ d(o)+Di{o)°^ an< ^ ^ an< ^ are un known functions. In the 
euclidean space, at the lowest order in perturbation theory one finds 

ZF-V)=0 Dr\x 2 )= 1 -^-^ 

07T X 

which is supposed to give the behaviour of the functions for x — > 0. A good 
parametrization for the long range region seems to be 

(3 D LR (x 2 ) = de _ff l*l, S = (1 ± 0.1) GeV , d = 0.073 GeV 4 , (30) 
f} D\ K {x 2 ) = d ie - Sllx \, Si = (1±0.1) GeV, di = 0.0254 GeV 4 , 

where-the value of gluonic condensate and lattice simulations have been com- 
bined llil. 



3.3 Dual QCD 

Dual QCDi is an effective theory described by a lagrangian density £ e g in 
which the fundamental variables are an octet of dual potential C 7 ' coupled to 
a classical quark source and to three octets of scalar fields Bi earring "mag- 
netic" color charge. The "monopole" fields Bi develop nonvanishing vacuum 
expectation values i?oi that give rise to a dual Meissner effect and provide a 
concrete realization of the Mandelstam-t'Hooft picture of confinement. We set 



where 



£ cff = 2Tr 



-~G^G^ + i(^B,) 2 



W(Bi), (32) 



X> M B< = 9 M Bi - i flM [C M ,Bi], (33) 
G>y = 9 M C W - ft,C M - i 9 M[C M , C v ] + G*,, (34) 
£gf is the gauge fixing term, 5 a/ = 2zc j s the magnetic coupling constant and 

= 9 ~k e * a 'i> a J da ^^r 5 ix-u{a,t)), (35) 

u' 1 = u M (i, a) being the surface swept by the Dirac string connecting the 
quark and the antiquark. The Higgs potential W(Bi) (which we do not give 
explicitly) has a minimum at non zero values Bqi of the form 

Bi = Bq\-j , B 2 = —B \ 5 , B 3 = B \ 2 ■ (36) 

In our present context the basic assumptioifl is that for large T 

W(T) ~ W cS . (37) 

To calculate explicitly W e ff even for restrict class of quantum fluctuation is 
not a trivial task. However, in the so called classical approximation this quan- 
tity is given by the expression exp{i J d 4 xCGp(G^ u )} evaluated for a classical 
solution of the field equations of the form 

Bi = BA7 B2 = — B\§ B3 = B A2 (38) 



where , B and B' are appropriate solutions of the equations 
^{d^-dyCy) = -d»G s ^-6g 2 C^B 2 

satisfying the asymptotic conditions — > 0, -B, B' =r* B$. Then Eqs. ^ can 
be solved numerically and analytically parametrizedQ. 

4 Semirelavistic Potential 

On the basis of general invariance principle the quark-antiquark potential up 
to order 1/m 2 can be written as 



^qq = V (r) + VvdO,Pi,p 2 ) + V r g D (r,pi,p2, 0-1,0-2) (40) 



with 



V V v(r(t)) = -J— ( Pl ■ p 2 H(r) + (\pi ■ p 2 - Pl " r 2 P2 " r ) I ; ( - ) 

mi 777,2 I \3 



Wcyl 



+ E ^ (^(r) + ( k 2 - Pj " r 2 Pi " r ) y.(r) } , (41) 
t J V3 r 2 7 J Woyl 



and 



Fsd = K4 + i )A[i;,,ri+v;(,)1 



^2 L i ' ff i " 'H;x W>(r) + 2V 1 (r)} 



imf Am\ ) r dr 

1 .Id 1 / (Ji • r (T2 • r 1 , 

(Li • 02 - L 2 • 01J -— K 2 (r) + 5 -o-i ■ o 2 K 3 (r) 



2mir7i2 r dr \m\mi \ r 2 3 

1 



1277717772 



a l -a 2 V i {r), (42) 



in terms of functions of r alone. 

ThepKajdous functions Vo(r), . . . V e (r) can be written in terms of momenta 
of fieldsO,El. Such expressions are very useful for lattice simulations!! 

Finally, one can work out a potential in explicit form for any of the models 
considered in Sec. 3. In Tab. 1 for the MAL model we report the complete 
expression, while for SVM and DQCD we report only the long range behaviour 
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Table 1: Complete MAL potential and long distance SVM and DQCD potentials 





MAL 


SVM 


DQCD 


V 


-l^ + ar + c 


cr 2 r+ iC*^ 1 ' -C 2 


err - 0.646 % 


/era s 


AV a 





Self-energy terms 






v{ 


a 






^crcts 


vl 


4 a B 
3 


r 






v 3 












Vi 


32-Ka a S' A (r) 










v b 


9 r 9 CT/ 




-0.097ar - 


0.226^/^07 


v c 


3 r 6 




-0.146ar - 


0.516^/<ra s 


v d 




l m . , V, 1^(1) ,10 2 2 B 2 


-0.118or + 


0.275v^al 






6 2 r 3 


-0.177or + 


0.258VctoI 



(we omit exponentially vanishing terms) . In terms of the parametrization of 
Eq. 30 the constants occurring in the second column are given bylj 

cr 2 = — C 2 = - (43) 

(1) 4di _ 3vrd 32d 

Complete expressions even for SVM and DQCD can be found in Eli andi. 
Notice that the short range behaviour must agree for the three potentials by 
construction. In the long range region, as one can see, there is essentially 
agreement among the leading terms. Some minor discrepancies in the DQCD 
case are possibly due to numerical inaccuracy. On the contrary the discrepan- 
cies in the subleading terms seem to be more significant. Notice however the 
complete similarities between the SVM and the DQCD columns for the first six 
lines of the table. Obviously the three models definitely differ at intermediate 
ranges. 

For what concerns the numerical values of the constants one may recall that 
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from heavy quarkonia fitting typical values are a ~ 0.18GeV 2 and a s ~ 0.37. 
On the other side if from lattice simulation we assume [cf. Eq.30] the marginal 
values 8 = 5! = 1.1, wc find a ~ 0.19GeV 2 , C 2 ~ 0.22GeV, cf } ~ 
0.08GeV, D 2 ~ 0.47, E 2 ~ 1.45GeV -1 . In this case the agreement be- 
tween SVM and DQCD is therefore remarkable for the first six lines even from 
the numerical point of view. Finally, notice that in the MAL case the spec- 
trum is very little sensible to the values of the constant C. In fact due to the 
combined form of Vq and Vn a variation of C can be reabsorbed in refinition 
of the quark masses m — ► m + This is not true in the SVM and DQCD 
cases due to the form of V\ and V 2 ■ 



5 Bethe-Salpeter kernel and effective mass operator 

Let us begin to notice that in the MAL case, neglecting the perimeter term 
and taking into account Eqs. [n] and 24, Eq. 19 can turn out to be of the 
general form of Eq. [l4| with only the first term. Additional terms would be 
necessary if we want to include higher order perturbation terms in Eq. |l9| or 
replace the simple MAL non perturbative part of the Wilson loop by more 
complicate expressions, like that of Eq. |37|. Of the terms in A the % = j terms 
correspond to self-energies of the quark or of the antiquark, the i =/= j terms 
to the interaction between the quark and the antiquark. 

Introducing Eqs. [l9[ 24 in 17,16 one finds explicitly in the momentum 
space (after factorizing the conservation delta (2-7r) 4 (5 4 (pi + p 2 — p[ — p 2 )) 



I(Pl, P2-, Pi, P 2 ) = Ipert + L 



167r-a s {iV(QK<72 



with 



16 



J dW Q r J(r, qi , q 2 ), 



2ar 



2 2 
9io920 

IqtI 



(arcsin 



QtI 



92o\/9io 



arcsm ■ 



|qri 

920 



2ar k <J^q 2oqil/ 



qio 

In such equations a t — — 

the free gluon propagator; furthermore we set qi 



(44) 



(45) 



02 giofe, 



4 — denotes the strong coupling constant and D pa (Q) 

^,92 



ind E 
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pi — p[ = p' 2 — p2 with qi = q 2 = q , = (S hk — f h r k )q k (CM frame has 
been assumed). 



5.1 Effective Mass Operator 

In terms of the relative and the total momentum the homogeneous equation 
corresponding to Eq. [l^ can be written 

* B (k) = -i J J^^ l \r 11 P B + k)HP( V2 PB - k)i(k,k';P B )$B(k'), (46) 

with r)j = ( mi " 3 m2 ) 5 Pb = (fiB,0), me being the mass of the bound state 
and $b an appropriate wave function. In the instantantancous approximation 
(consisting in setting Hj(p) = - » , and k Q = k' = r\ 2 Wl+ 2 Wl - r/i ™ 2 ^™ 2 in 



I(k,k';P) with = ^/m 2 + k 2 , = ^/m 2 + k') the residual variables /c 

and fcg can be integrated explicitly and one is left with the eigenvalue equation 
for an effective mass squared operator M 2 — Mq + U with Mq — wi +w 2 and 



WHO = ^R^^M%^. (47) 



(27r) 3 V 2wiw 2 v ' 'y 2w'iw' 2 

From this one can also derive the more conventional operator M = Mo + V, 
where 

(k\V\k>) == 1 r J ==i inst (k,k>) + ... (48) 

having kept only first order terms in U and neglected kinematical factors equal 
to 1 on the energy shell (wi + w 2 = w[ + w 2 ). 
From Eqs. 45 one obtains explicitly 



<k|[/|k') = 



(wi +w 2 )(w[ + w' 2 ) ( 4a s r 1 r 2 (Q • q 



WiW 2 w'iW 2 \ 3 7T 2 



,2 



Q -L<Zl0920 + q- , Q2 



i 1 

k X k' • (ai + (T 2 ) + -TT^lqmiai ■ Q) - q w (a 2 ■ Q)] 



2Q 2 v " 2Q 2 

\°i • °* + • °* - (Q '^ ( 2 Q " J2) ) + • Q)(a 2 • Q) (49) 



(20 



^3 | d 3 re^ r J inst (r,q, 9lo ,g 2 o) 
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with 



J inst (r, q, q 10 , q 20 ) = ■ [q% J q 2 - q 2 + q 2 J q 20 - q 2 ,] + 

9l0 + 920 v v 

910920/ • I<3t| . . |qTL 
+ (arcsm h arcsm )J 

|Qt| 9io 920 

a ! 920 / , . / s s . 9i0/ / 
i / 2 2 : ( r x q • ^1 + i9i (r • ai)) + (r x q • a 2 - i92o(r 

r V 9io _ qT V 920 - qT 



Here denote the usual Dirac matrices 7^7^, ff* the 4x4 Pauli matrices 
and obviously q = "t , Q = k k', 9jo = 3 2 J . Notice 





a k 



that, due to the terms in ctg, such {/ is hermitian only with reference to the 
undefined metric operator 7i72- 



Due to Eq. 48 the potential V can be obtained from U as given by Eqs. 



simply by the kinematical replacement 



(wi+w 2 )(w' 1 +w' 2 ) 1_ 



w 1 w 2 w' 1 w 2 2y / w 1 w 2 w' 1 w' 2 ' 

Notice that, if in the resulting expression we perform an 1/m expansion and 
make an appropriate Foldy-Wouthuysen transformation to eliminate the terms 
, we reobtainB the MAL potential as given in Tab. 1. 
On the contrary, if we keep only the long range terms in V , neglect the 
spin dependent part and set for simplicity mi = m 2 = m, we can write in an 
operatorial form 



m a k 



M = 2^/m 2 + q 2 




(Ql = q _ qT) in which the ordering implied by the definition referring to Eq. 50 
has to be understood. Such expression is identical to the hamiltonian of the 
Rclativistic Flux Tube Model 13 up to the order era 2 . Two different limits 
of |5l] are of interest. The first one is the limit of small angular momentum, 
i.e. small transversal momentum (q^ = -^j-). In this case we have (strictly for 
s waves) 

M = yjm? + q 2 + or (52) 

This result justifies the use that of Eq. [52| has been done in the study of 
light meson spectrum. The second limit is for large angular momentum or 
transversal momentum (negligible q 2 ) and is 



M = 2V^T^ + ? ( V^I+Z arcsin \*L= + -f^—) (53) 
2 |q| yjm 2 + q 2 Vm 2 + q 2 



13 



or under the additional assumption of negligible m 

M = 2|q| + ^ar (54) 

Eq. ^4] is very important for an understanding of the Regge trajectory prop- 
erties. In fact it is well known that such hamiltonian produces asympotically 
straight trajectories with 

mf 8^<rl = 2ircrl (55) 

mi being obviously the mass of a given radial quantum number bound state 
as a function of the angular momentum I. Notice that Eq. |55| corresponds to 
the slope a' = -fi? — > 7^ which is identical to that of the Nambu-Goto string 
model. Notice that for a = 0.18 we find 0/ = 0.88 in very good agreement 
with the the experimental slope of the p trajectory. Had we used the naive 
Eq. [5^, even for large I, we would have find a' = l/8a p. 0-69. As mentioned 
such results are in perfect agreement with those of Refsll3,ll3. 

Finally, if we keep all terms in Eq. but assume ni2 3> ni\, we can discuss 
the light-heavy quark symmetry much on the lines followed by KaidalovliS. 
E. g. let us consider the quadratic hyperfine separation between triplet and 
singlet 5m 2 = m\ — m| in the (uQ) states. Empirically this quantity is nearly 
independent of Q (~ 0.55 MeV 2 ). In fact the hyperfine splitting term in 
U (cf. [19]) depends on the quark masses only through the kinematical factor 
occurring in it. Now for m 2 » ni\ such factor reduces to . — \ , and any 
dependence on nii disappears. 
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